Even-odd entanglement in boson and spin systems 
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We examine the entanglement entropy of the even half of a translationally invariant finite chain 
or lattice in its ground state. This entropy measures the entanglement between the even and 
odd halves (each forming a "comb" of n/2 sites) and can be expected to be extensive for short 
range couplings away from criticality. We first consider bosonic systems with quadratic couplings, 
where analytic expressions for arbitrary dimensions can be provided. The bosonic treatment is then 
applied to finite spin chains and arrays by means of the random phase approximation. Results for 
first neighbor anisotropic XY couplings indicate that while at strong magnetic fields this entropy is 
strictly extensive, at weak fields important deviations arise, stemming from parity-breaking efi'ects 
and the presence of a factorizing field (in which vicinity it becomes size-independent and identical 
to the entropy of a contiguous half). Exact numerical results for small spin s chains are shown to 
be in agreement with the bosonic RPA prediction. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 

The entanglement properties of many-body systems 
are of great interest for both quantum information theory 
[l[ and condensed matter physics [343| ■ Their knowledge 
enables, on the one hand, to assess the potential of a 
given many-body system for quantum information pro- 
cessing tasks such as quantum teleportation ^ and quan- 
tum computation 1, 6, 7|]. On the other hand, it provides 
a deep understanding of quantum correlations and their 
relation with criticality [2-0, S 0|- In non-critical sys- 
tems with short range couplings, i.e., local couplings in 
boson or spin lattices, ground state entanglement is be- 
lieved to satisfy a general area law by which the entropy 
of the reduced state of a given region, which measures its 
entanglement with the rest of the system, scales as the 
area of its boundary as the system size increases 0, fioj . 
This behavior is quite different from that of standard 
thermodynamic entropy which scales as the volume. In 
one dimensional systems this statement has been quite 
generally and rigorously proved 0, [ll| and simply means 
that the entropy of a contiguous section saturates, i.e., 
approaches a size independent constant, as the size in- 
creases. Violation of this scaling is therefore an indica- 
tion of criticality |^, . The exact expression of the 
entropy of a contiguous block in a one-dimensional XY 
spin 1/2 chain in the thermodynamic limit has been ob- 
tained and confirms the previous behavior. 

The conventional area law holds for contiguous sub- 
systems. For non-contiguous regions it actually implies 
that the entropy is proportional to the number of cou- 
plings broken by the partition. For instance, for comb- 
like regions like the subset of all even sites in a chain, the 
entropy should scale as the total number n of sites for 
first neighbor or short range couplings. This was in fact 
verified in HT] for the harmonic cyclic chain, where the 
corresponding logarithmic negativity was calculated, and 
also verified numerically in [16[ for some spin arrays and 
a 1-d half-filled Hubbard model, where the even entan- 



glement entropy was computed. An exact treatment of 
general comb entropies for a large one-dimensional criti- 
cal XX spin 1/2 chain with first neighbor couplings was 
given in [IT] , showing that they are indeed proportional 
to the size L plus a logarithmic correction. 



The aim of this work is to analyze in detail the en- 
tanglement entropy of all even sites in finite boson and 
spin arrays, both in one dimension as well as in general 
d-dimensions. Such bipartition can be normally expected 
to be the maximally entangled bipartition at least for uni- 
form nearest neighbor couplings, as it will there break all 
coupling links. We first analyze the bosonic case with 
general quadratic couplings, where a fully analytic treat- 
ment of this entropy is shown to be feasible and allows 
to derive simple general expressions in the weak coupling 
limit. Comparison with single site and block entropies 
is also made. The bosonic treatment is then applied to 
finite spin s arrays with anisotropic ferromagnetic-type 
XY couplings in a uniform transverse field through the 
RPA approach fl^. This allows to predict in a simple 
way the main properties of the total even entropy in these 
systems. Comparison with exact numerical results indi- 
cate that the RPA prediction, while qualitatively correct, 
is also quite accurate outside the critical region already 
for low spin s > 2, representing the high spin limit. Re- 
sults corroborate that for strong fields, the total even 
entropy in these systems is extensive, i.e., directly pro- 
portional to the total number n of sites. However, for 
low fields B < Be, this entropy has an additive con- 
stant, which arises in the RPA from parity restoration 
p^ . Moreover, in the immediate vicinity of the factoriz- 
ing field Bs < Be 19,-22] , extensivity is fully lost and the 
total even entropy reduces to this constant, which is the 
same as that for the block entropy and is exactly eval- 
uated. The exact bosonic treatment is described in sec. 
im whereas its application to spin systems is discussed in 
sec, mil Conclusions are finally drawn in IIVI 
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II. ENTANGLEMENT ENTROPY IN BOSONIC 
SYSTEMS 

We start by considering a system of n bosonic modes 
defined by boson creation operators bj {[bi,b'j] = 6ij), 
interacting through a general quadratic coupling. The 
Hamiltonian can be written as 

H = - ^)(blb, + ^S,,) ~ ^AT^blb] + A7.bjl>\) 

-A- 



where Z''^ = {b\b), Ay — XiSij and the 2n x 2n matrix 
H is hermitian. The system is assumed stable, such that 
the matrix "H is positive definite. We may then also write 
^ in the standard diagonal form 



which has eigenvalues and —1 — /^), and given by 



S{pa) = -TrpAlnpA ^}_^hi.f^) , (6) 



k=l 



M/) = -/ln/+(l + /)ln(l + /) 



(7) 



For instance, the entanglement of a single mode i with 
the rest of the system is just 



S{p^) = h{f,) , /. = ^{F+ + ^r-\Fr\^ - i , (8) 

where fi, the symplectic eigenvalue of the single mode 
contraction matrix 2?^, represents the deviation from 
minimum uncertainty of the mode: {F^ + i)^ — 1^7 P ~ 



{qf)o'{pl)o' - [Re((g.p.)o')]' > for q, 

\/2'i ■ 



bi+b 
~72 



Pi = 



H 



E 

k 



Mb'lb', + i) 



(2) 



where cok are the symplectic eigenvalues of i.e., the 
positive eigenvalues of the matrix MH, with M = ), 
which come in pairs of opposite sign and are all real non- 
zero when y. is positive definite [11], and b'l are the 
normal boson operators determined by the diagonaliz- 
ing Bogoliubov transformation Z — WZ' satisfying 
W^'MW = M and {W'''HW)kk' = t^kSkk'- The ground 
state is the vacuum |0') of the operators 6J, and is non- 
degenerate. 

Ground state entanglement properties can be eval- 
uated through the general Gaussian state formalism 
prl . [23 , [25I I , which we here recast in terms of the con- 
traction matrix [isl l23t 



A. Finite translationally invariant systems 



Let us now associate each bosonic mode with a given 
site in a cyclic chain and consider a translationally in- 
variant system of n sites, such that Ai = A and = 
A^{i - j), with A±(-;) A^{n~ I). We first con- 
sider for simplicity the one-dimensional case. Through 
a discrete Fourier transform &| — X]fc=o '^'^'^'^'^"^L 
we can diagonalize analytically and obtain an explicit 
expression for the contractions F^. We will assume 
A^(/) = A^{—1) V I, in which case the energies uik in 
dl]) adopt the simple form [3| 



^k = x/{x^A+y-{A-y, 



(9) 



where Aj, are the Fourier transforms of the couplings: 



V = (ZZ^)o' -M=W 



F^ 
F- 



F-_ 





1 ^ 

F+ = {b]b,)o' 

T = {b,b,)o, = {bib])*, 



(3) 



This hermitian matrix determines, through application of 
Wick's theorem [2^ , the average of any many-body oper- 
ator. In particular, the reduced state pA = Tr^|0')(0'| of 
a subsystem A of ua modes (A denoting the complemen- 
tary subsystem and Tr^ the partial trace) is fully deter- 
mined by the corresponding sub-matrix Da = {ZaZ\) — 
Ma (Eq- Q with i,j e A) and can be written as [ll| 



PA 



exp[-^z\HAZA]/Tre^p[~^^z\HAZA] , (5) 



where Ha = Ma^t^[I + Ma'Pa^]- Eq. ^ represents a 
thermal-like state of suitable ua independent modes de- 
termined by the effective Hamiltonian 'Ha- The entan- 
glement entropy of the {A, A) partition, S{pa) = S{pa), 
is then determined by the symplectic eigenvalues of 
Da (i.e., the positive eigenvalues of the matrix VaMa, 



A^ ^ ^ei27rfc;/"A±(0. 



(10) 



(4) The contractions F^ depend just on the separation 



\i — j\ and are given by 

n — 1 



(11) 



fc=0 



n = {bibkW = i = {bkb-k),' = ^2) 

The symplectic eigenvalues of the full contraction matrix 

V fc. 



p-\2 _ 1 



Q are of course fk = ^ {\ + f^f - 

In the weak coupling limit |A^| ^ A V fc, become 
small and up to lowest non-zero order we obtain 



fk 



a: 



_ (A^ 
2A ' -"^ ^ 4A2 ^ ' ' 



' J k 



if, 



-^2 



(13) 



which leads to 



Fr.^.F,^.^lL^:iDS^. (14) 
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FIG. 1. (Color online) Even-odd partitions of one- and two- 
dimensional arrays 



At this order just sites linked by A^(Z) or its convolution 
are correlated. The eigenvalues of subsystem contrac- 
tion matrices will depend up to lowest non-zero order on 
F+ and (Ff )2, being then ©(A?. /A^) for A-{1) cx A_. 
We can then use in ([6]) the approximation 

/i(/)«-/(ln/-l) + 0(/2), (15) 

such that S{pa) — In -jr)- 

On the other hand, it is seen from Eq. ^ that the 
present system is stable provided A > A^ -I- |A^| V k. 
For attractive couplings A+(Z) > V I, with all A~(Z) 
of the same sign, the strongest condition is obtained for 
fc = 0, so that stability occurs for 

A > A, = A+ + |A„ I = A+{1) + \A-{1)\ . (16) 

I 

For A — > Ac, Wo — (while all other Wfc remain finite in 
a finite system), implying a divergence of Jq (Eq. ([12]) '): 



\fo 



8(A-Ac)' 



/o+« 1/(71-1/2, (17) 



plus terms 0{X/Xc — 1)^/^. This entails in turn a di- 
vergence Jq oc (A/Ac — 1)^^/^ of the largest eigenvalue 
of a subsystem contraction matrix V^, with S{pa) ~ 
In + 1 !v — i ln(A/Ac — 1) plus constant terms. 
For example, the single site entropy ([5]) becomes 



S{p^)^Hf), f^.^(l+F+r-iF,-r-^, (18) 



with F(± = i Efe fk (Eq- dl))- For weak coupling. 



-^2 



4A2 



(19) 



which involves just the couplings A (I) connecting the 
site with the rest of the system. On the other hand, for 

A^Ac, /cx v^(x(A/Ac-l)-i/4. 



B. Even-odd entanglement entropy 

We now evaluate the entropy of the reduced state of 
all even sites, S{pe) = S{po), which measures their en- 
tanglement with the complementary set of odd sites (Fig. 



[l]left). We will assume n even, such that the even sub- 
system, defined by (—1)* — +1, is again translationally 
invariant. The ensuing contraction matrix Ve can be ob- 
tained by removing contractions between even and odd 
sites in the full matrix ([4]) and extracting then the even 
part. This leads to elements 



f;± = i^;±(i-Fe-(-^")), 

whose Fourier transforms are, using Eq. (Ilip . 



(20) 



~fk=Wt^ftn/2)- (21) 

The final symplectic eigenvalues of 25b then become 



for fc = 0, . . . , n/2 — 1. We then obtain 



n/2-l n-1 

s{pe) = '^(^) = 5 E '^(/fe) 

fe=0 fc=0 



(22) 



(23) 



Whenever can be approximated by a smooth function 
/(fc) of A: = fc/n, we may replace ((23l) by the integral 



S{pe) 



h[j{k)]dk. 



(24) 



In these cases, we may then expect S{pe) extensive, i.e., 
proportional to the number n/2 of even sites. Let us re- 
mark, however, that this is not always the case: In a com- 
pletely and uniformly connected system like the Lipkin 
model [13, m, 113] , the contraction matrix will have a sin- 
ofe non-zero symplectic eigenvalue for any subsystem 
p^ . including the whole even set, and S{pe) = h{fn/2) is 
no longer proportional to n. A similar lack of extensivity 
holds in a finite system in the vicinity of the instability 
(A — >■ Ac, see below). 

For weak coupling, Eqs. (fT3|) . (fT4|) and (f2T|) lead to 



fk 



(A, 



^'k+n/2f (E;odde^'^'^'/"A-(/))^ 



16A2 



4A2 



(25) 



which involves again just the couplings A (Z) connecting 
the even and odd subsystems. On the other hand, for 

A ^ Ac (Eq. dm), /o « i^(l + 2/+2)|/o-|-2/-/2/o" 

diverges as (A/Ac ~ 1)^^^"' whereas all other f^^ remain 
finite, and extensivity is lost. 



C. First neighbor coupling 

Let us now examine in detail the first neighbor case 
A±(0 = \A^{5ii + 6i-i), where Eq. ^ becomes 



A^ = A=^ cos(27rfc/n) 



(26) 
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The exact S{pe) can be obtained from Eqs. ((2T|) -(l23l 
In the weak coupHng hmit, Eqs. ([TS]) and lead to 



(A-)^ 
8A2 



/fe«%^~2/cos2(2^A;/7i). 



(27) 
(28) 



Using Eqs. (|15I) - ([M|) . the single site and the total even 
entropies can then be expressed just in terms of /: 



5(p,)«-/(ln/-l) 



(29) 



S{pe) « -nf /cos2(27rfc){ln[2/cos^(27rA:)] - 1} dk 



--/(ln/-ln2). 



(30) 



Hence, in this limit S{pe) is extensive^ becoming n/2 
times the single site entropy (1291) minus a 0{nf) correc- 
tion accounting for the interaction between even sites: 



(31) 



The last term represents the even mutual entropy 
^S{pi) — S{pe), which is always a positive quantity and 
becomes here also extensive in this limit. 

In contrast, the block entropy S{pl), where pL denotes 
a contiguous block of L < n spins, rapidly saturates as 
L increases In the weak coupling limit, it is ver- 



ified that the ensuing contraction matrix Vl possesses, 
up to lowest non-zero order, just two positive non-zero 

i/ for any L >2, such that 



symplectic eigenvalues 



S{pl) ~ -/(In f/2 - 1) « S{p,) + / In 2 , (32) 

for 2 < L < n — 2, i.e., it saturates already for L = 2. 
Hence, in this limit. 



S{PE)^^SipL)-^f. 



(33) 



Assuming A+ > (if A+ < we can change its sign by 
a local change bi — >■ —bi at odd sites) the present system 
is stable for A > Ac = A+ + |A_| (Eq. (US])). For A ^ A^, 
Wo — and all previous entropies diverge. In particular, 
Eq. dUl) leads to 




|A-|Ae 



23A+(A- Ac) 



-1] 



being then verified that S{pe) ~ —j ln(A/Ac — 1) plus a 
constant term up to leading order. Hence, in this limit 
S{pE)/Sip,) ^ 1. 

As illustration, the left panels in Fig.[5]depict the single 
site, block and even-odd entanglement entropies for a ring 
of n = 36 sites with A" = A+/3, where Ac = 4A+/3. 



D. Even-Odd entropy in d-dimensions 

The whole previous treatment can be directly extended 
to a translationally invariant cyclic array in d dimen- 
sions. We should just replace l,k,n by vectors I = 
{h,.. Id), k = (fci, . . . , fed) and n = (ni, . . . , n^), with 
Zi, fci = 0, . . . , TT-i — 1. We will assume couplings satisfying 
A±. = A±(2 - j), with A±(-0 = A±(n - I) ^ A^{1). 



The same previous expressions (|9))- (|T2t then hold, with 
A± =^e*2-^-'A±(Z), (34) 



(35) 



where k = [ki/ni, . . . ,kd/nd) and n = JliLi '^i is the 
total number of sites. Eqs. (|13p - p^ remain unchanged 
with i,k,l ^ i, k, I. 

The subsystem of all even sites, like that formed by 
the blue sites in Fig. [T] right, is defined by 

Its contraction matrix will then be the even block of 

where 1 = (1,...,1). Assuming Ui even V i, its Fourier 
transform is then given again by 



(37) 



where ki + ni/2 — >■ ki~ni/2 if fc^ > ni/2. The symplectic 
eigenvalues of "De are then given again by Eq. (|22p with 
k ^ k, and the even-odd entanglement entropy reads 

SiPE) = \ E '*(/'^) ^ f / h[~f(k)]d''k , (38) 

k 

where ki = 0, . . . , — 1 in the sum and the integral is 
restricted to the unit cube < fc^ < 1 and valid if fk is a 
smooth function f{k) of k. 

In the case of first neighbor couplings 

A±(Z)^i^Af(5^e.+5L-eJ, 



where e.^ = (0, . . . 1^, . . . 0), Eq. ([Ml) leads to 

d 

A± = E A± cos(27rfc,/n,) . 



(39) 



with A^^^y2 = ^AJ. In the weak coupling limit we 
then obtain 



8A2 ' 



i=l 



(40) 



/.«u(fe)/, u(fe) = 2(E^cos^)^ (41) 
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FIG. 2. (Color online) Top left: Entanglement entropies in the ground state of a one-dimensional bosonic chain of n = 36 
sites described by the Hamiltonian ([1} with first neighbor cyclic couplings and A" = S{pi), S{pl) and S{pe) denote, 

respectively, the entropy of the reduced state of a single site, a block of L — n/2 contiguous sites and the set all even sites. 
Bottom left: The ratios S{pe) / S{pi) and S{pl) / S{pi). Dotted lines depict the ratios determined by the asymptotic expressions 
(|29[) - (|32|) . For large A these ratios are then close to n/2 and 1 respectively, while for A — > Ac they all approach 1. The right 
panels depict the same quantities for a two-dimensional square array of n = 6 x 6 sites with isotropic cyclic couplings and the 
same ratio A~/A^. S{pl) denotes the entropy of a contiguous half of 6 x 3 sites. Dotted lines correspond to the expressions 
At fixed A/Ac, S{pe) and S[pi) are now roughly half the value of the left panel (Eq. (gS))), while S{pL)/S{pi) is 
proportional to ^/n (Eq. (|46|) ). 



Hence, the single site entropy is agam S{pi) 
1) while Eq. yields 



--/(ln/-l + a) 



^./(In./^ decreases: 



S{pe) 



where a is a geometric entropy factor: 
u{k)lnu{k)d'^k, 



(42) 
(43) 

(44) 



(w(fc) > 0, J u{k)d'^k — 1). In the isotropic case A,^ = 
A~ V i, we have a = ad, with ai = 1 — ln2 w 0.307 
(Eq. [nj), a2 — 2q!i « 0.614 and 0:3 ~ 0.636, approaching 
« In 2 for large d. 

At fixed A, and for Af = A±, / = {A-)^d/{8\^) 
and hence both S{pi) and S(pe) increase as d increases, 
reflecting the larger number of links. However, and as- 
suming again A+ > 0, Ac = d(A+ + |A~|) also increases, 
entailing that at fixed A/Ac, / (and so S{pi) and S{pe)) 



f 



[A-/(A^ 



8d(A/Ac 



(45) 



For example, the right panels in Fig. [5] depict S{pe) and 
S{pi) in an isotropic square lattice of 6 x 6 sites, with 
the same previous ratio A~/A+ = 1/3. At fixed A/Ac, 
their values are verified to be roughly half that of the 
similar one-dimensional case (Eq. (|1S])). Their ratio is 
also slightly smaller due to the increase in the parameter 
a in P5|) . On the other hand, for A — Ac there is again 
a single vanishing energy ujo, so that all entropies behave 
as — iln(A/Ac — 1) up to leading order, with all ratios 
approaching 1. 

We also depict there the entropy S{pl) of a contiguous 
half-size block (n^ x ny/2 = 6x3 sites), which is now 
proportional to its boundary 2nx. For A ^ Ac, it is ver- 
ified that the number of non-zero positive eigenvalues of 
the corresponding contraction matrix is just the num- 
ber of couplings "broken" by the partition {2nx), being 
all approximately equal to //4 up to leading non-zero 



6 



order. We then obtain 

77 71 

S{pl) « - Y/(ln//4 - 1) « -f[S{pi) + 2/ In 2] , (46) 

whence S{pl)/ S{pi) oc nx/2 in this hmit, as verified in 
the right panels of Fig. [51 



III. APPLICATION TO SPIN SYSTEMS 

The previous bosonic formalism can be directly applied 
to interacting spin s systems in an external magnetic field 
through the RPA approximation [isj . Denoting with s^^ 
the dimensionless spins at site i, we will consider a 

cyclic translationally invariant finite array which can be 
described by an XY Hamiltonian of the form 

H = B^s^z- ^ Y^i-J^jSi^Sj.^ + JfjS^ySjy) (47a) 



where Sj± = Sjx ± isjy, J^J — Jfj.ii ^ j) and 



(48) 



We note that x, y, z may in principle also denote local in- 
trinsic axes at each site, in which case the field is assumed 
to be directed along the local z axis. The scaling of 
the couplings ensures a spin- independent mean field and 
effective RPA boson Hamiltonian (see below). 

Normal RPA. For sufficiently strong field B, the low- 
est mean field state (i.e., the separable state with lowest 
energy) is the aligned state |0) = |0i) ® . . .® |0„), where 
|0i) denotes the local state with maximum spin along the 
— z axis {siz\Qi) = —s\Qi)). In such a case, RPA implies 
the approximate bosonization [l8l | 



2sb, 



b]b,. 



(49) 



which is similar to the Holstein-Primakoff bosoniztion 
[23I [2^ and leads to the quadratic boson Hamiltonian 
([T]) with the parameters (j48|) and X = B. We may then 
directly apply all previous expressions. 

The bosonic RPA scheme becomes exact for strong 
fields \B\ ^ Be for any size n, spin s, geometry or in- 
teraction range, since for weak coupling it corresponds 
to the exact first order perturbative expansion of the 
ground state wave function As a check, in the case 
of the spin 1/2 one-dimensional chain with first neighbor 
XY coupling, an analytic expression of the block entropy 
in the limit n — > cxi has been obtained in [13l4l5i . For 
X = B > A+ , it is given in present notation by [l3| 



,2x2£(a)JV)i 



(50) 



where a = A'/Va^ -I- A"^ - A+^, a' = Vl - and 
I (a) = dx/ ^/(l — a;^)(l — a^a;^) is the elliptic integral 



of the first kind. An expansion of ([50)) for A ^ A^ leads 
exactly to present Eq. ([5^ . with / given by ([?7)) . We can 
then expect the asymptotic expressions ([30)) and ([43]) for 
S{pe) to be exact in this limit also in spin systems. 

Parity breaking RPA. Considering now the anisotropic 
ferromagnetic- type case |Jy(OI < Jx{l) V / in (|47ap . the 
previous normal RPA scheme will hold, according to Eq. 
([THj), for B > Be = = Jx{l), i-e., when the corre- 
sponding boson system is stable. 

For \B\ < Be, the normal RPA becomes unstable {luq 
becomes imaginary). The lowest mean field state corre- 
sponds here to degenerate states | ± 8) fully aligned along 
an axis z' forming an angle ±6 with the z axis in the x, z 



plane: |e) = \0i) 



{On), with \0,) = exp[-i9siy]\0i). 



We are assuming here an anisotropic XY coupling such 
that H commutes with the Sz parity P^ = e*'^^^* si^+ns) ^ 
but not with an arbitrary rotation around the z axis (as 
in the XX case). Such states break then parity sym- 
metry, satisfying Pz\'c>) — \ — Q). The angle 9 is to be 
determined from llSlI 



cosd^B/Bc, Be 



I 



MO- 



(51) 



For |B| < Be, the bosonization (|49p is then to be applied 
in the RPA to the rotated spin operators Siz' — Siz cos 9+ 



. sin 6*, Si±' = Six'±isiy', with Si^' = s^^r cos 



; sm ( 



and Siyi — Siy. This leads again to a stable Hamiltonian 
of the form ([T]) with aS] 



X^Be, A^{1) = ^[r{l)cos^9±jy{l)]. 



(52) 



For \B\ < Be we should also take into account the im- 
portant effects from parity restoration for a proper RPA 
estimation of entanglement entropies [3]. The exact 
ground state in a finite array will have a definite parity 
Pz outside crossing points [2l|, implying that the actual 
RPA ground state should be taken as a definite parity 
superposition of the RPA spin states constructed around 
I ± 6) [l8|. This leads to reduced RPA spin densities 
of the form pA ~ ^[pa{(^) + PA{—d)] if the complemen- 
tary overlap — (^©aI^a) can be neglected. If the 
subsystem overlap Oa — (— ©aI©^) = cos"^"^^ 9 is also 
negligible, such that pa{9)pa{—9) « 0, then [ll] 



S{pa) ~ S{pa{9)) + 5 , 



(53) 



where 5 = In 2. The final effect is then the addition 
of a constant shift to the bosonic subsystem entropy for 
\B\ < Be. This is applicable to both S{pe) and S{pl) if 
9, n and the block size L are not too small. 

For first neighbor couplings with anisotropy x — 
Jy/Jx G (0,1) (if X > 1 we just redefine the x,y axes) 
as well as for arbitrary range couplings with a common 
anisotropy x — Jy{0/Jx{l) G (0,1), another fundamen- 
tal feature for \B\ < Be is the existence of a transverse 
factorizing field Bs — Be^/x where the mean field states 
I ± 9) become exact ground states [l9l - l22j | . As seen from 
([5^ . at this field A-{1) = V /, so that the RPA vac- 
uum remains the same as the mean field vacuum llSl and 
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FIG. 3. (Color online) Top: Exact entanglement entropy of all even sites (left) and of a contiguous block of n/2 sites (right) 
in the ground state of a one dimensional cyclic chain of n = 8 spins with anisotropic XY first neighbor couplings {Jy/Jx ~ ^) 
and spin s — 1/2, 1 and 2, as a function of the transverse magnetic field. The dotted line depicts the bosonic RPA result, with 
Be — Jx the mean field critical field. We have used base 2 logarithm in the entropy, such that all entropies approach 1 at the 
factorizing field Bs ~ 0.71Bc- Bottom: Left: The corresponding ratio S{pe) / S{pl)- Right: The entanglement entropy of all 
even sites in a rectangular lattice of 4 x 2 spins. Remaining details as in the top panels. 



all contractions vanish, implying S{pa{6)) — 0. All at Bg [Ic 



RPA entropies at Bg reduce then to the correction term 
5 arising from parity restoration . 

This is essentially also the exact result at Bg'- The 
transverse factorizing field corresponds to the last ground 
state parity transition as B increases from i21j and the 
ground state side-limits for B ^ Bg are actually the def- 
inite parity combinations of the mean field states | ± 8) 
[21, 28]. These definite parity states have Schmidt num- 
ber 2 for any bipartition, implying that the side-limits of 
the exact entropy of the reduced state of any subsystem 
at Bg do not approach but rather the values [l^ 



S(pJ(B.)) 



2(1 ±OaO, 



(54) 

where -I- (— ) corresponds to positive (negative) parity, i.e. 
the right (left) side limit at Bg [28]. Eq. dSH) is valid _for 
any size or spin. For small complementary overlap O^, 
~ 5(1 + a) and both side limits coincide, while if 
Oa is also small, « 1/2 and Eq. ((M)) reduces to -I- In 2. 
This is also in agreement with the exact limit of the block 
entropy of the large one dimensional s = 1/2 XY chain 



Illustrative exact results for the even-odd entanglement 
entropy in a finite linear cyclic spin s chain with first 
neighbor couplings are plotted in the top left panel of Fig. 
[3] for spins s — 1/2, 1 and 2, together with the bosonic 
RPA estimation. The exact definite parity ground state 
was employed in all cases. We also depict for comparison 
the entropy of a contiguous half (top right), and the ratio 
S{pe)/S{pl) (bottom left). The anisotropy of the cou- 
pling is the same as in Fig. [U (A_/A+ = 1/3). The RPA 
result (independent of s for the scaling used in (|T7)) ). rep- 
resents the large spin limit but is already quite close to 
the exact results for s = 2 except in the vicinity of Be, 
where the exact entropies remain of course finite in a fi- 
nite chain. The ratio S{pe) / S{pl) is nonetheless quite 
accurately reproduced and shows the extensive character 
of S{pe) for B > Be, in agreement with ([55)1 . where the 
entropies for all spin values rapidly approach the RPA 
result and become spin independent. For \B\ < Be the 
shift 6 in (|53p (^ = -1-1 in Figs [SHU since base 2 logarithm 
was employed) is essential for the agreement and explains 
the lack of direct extensivity in this region. The collapse 
of all entropies to the value S at Bg is also verified, and 
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1 2 

B/Be 




B/Be 

FIG. 4. (Color online) Top: Exact entanglement entropy of 
all n/2 even sites in the ground state of a spin 1/2 cyclic chain 
for different values of n. Couplings are the same as in Fig. [S] 
Bottom: The intensive entropy S(pE)/nE {jie ~ n/2). All 
curves coalesce for B > A+. The inset depicts the intensive 
shifted entropy {S(pe) — 5)/nE, where 5 = for B > Bs and 
5 = 1 for B < Bs, which makes curves for n > 8 coalesce also 
for B < Bs. 

for s = 1/2 even the small discontinuity at Bs predicted 
by Eq. ([5i)) can be appreciated (together with the other 
parity transitions for B < Bg). The bottom right panel 
depicts S'(p£;) in a 4 X 2 square lattice with identical 
couplings in both directions and the same ratio A_/A+, 



where a similar behavior is obtained. Exact results for 
.3 = 2 are now even closer to the RPA prediction, indi- 
cating that the accuracy of the latter tends to improve, 
for stable mean fields, as the connectivity increases poj. 

Exact results for a spin 1/2 chain for different sizes 
are depicted in Fig. Even though RPA is not accurate 
for such low spin with a first neighbor coupling [2^, the 
exact results are again in qualitative agreement with its 
predictions away from the critical region: Direct exten- 
sivity S{pe) oc n is verified for strong fields B > A-(_ 
(bottom panel), whereas for B < Bs it holds for the 
shifted entanglement entropy S{pe) — S, as seen in the 
inset. Complete lack of extensivity takes place at the fac- 
torizing field Bg , where the discontinuity implied by ([5^ 
is appreciable for ti = 8 and becomes quite noticeable for 
n = 4. 

IV. CONCLUSIONS 

We have shown that the total even-odd entangle- 
ment entropy displays a strict extensive behavior in both 
bosonic and spin chains or lattices for weak first neigh- 
bor couplings (i.e., strong fields in a spin chain), pro- 
viding explicit asymptotic expressions for the general d 
dimensional case. Extensivity of the associated mutual 
information is also implied by these expressions. Devi- 
ations from this behavior, however, were shown to arise 
for stronger couplings, i.e., proximity to the instability in 
the finite bosonic case or low fields in the spin case. In 
the latter, a constant shift is essential to understand the 
exact results for \B\ < Be, which has an evident meaning 
as a symmetry restoration effect in the RPA. Besides, full 
loss of extensivity occurs in the vicinity of the factorizing 
field. Present results confirm the validity of the RPA ap- 
proach (with inclusion of symmetry-restoration effects) 
for obtaining a simple direct understanding of the main 
aspects of ground state entanglement in spin chains, at 
least in those regions where a well defined mean field 
minimum exists. 

The authors acknowledge support from CIC (RR) and 
CONICET (NC,JMM) of Argentina. 
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